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Abstract 

We suggest a perturbative approach for generic choices for the universe equation of state and 



o 



O 

' introduce a novel framework for studying mass varying neutrinos (MaVaN's) coupled to the dark 

c/3 '. 

sector. For concreteness, we examine the coupling between neutrinos and the underlying scalar 
! field associated with the generalized Chaplygin gas (GCG), a unification model for dark energy and 

> 

' dark matter. It is shown that the application of a perturbative approach to MaVaN mechanisms 

m ■ 

' translates into a constraint on the coefficient of a linear perturbation, which depends on the ratio 

, between a neutrino energy dependent term and scalar field potential terms. We quantify the effects 

' on the MaVaN sector by considering neutrino masses generated by the seesaw mechanism. After 

o ■ 

setting the GCG parameters in agreement with general cosmological constraints, we find that the 
^ ■ squared speed of sound in the neutrino-scalar GCG fluid is naturally positive. In this scenario, the 

H : 

5^ , model stability depends on previously set up parameters associated with the equation of state of 

the universe. Our results suggest that the GCG is a particularly suitable candidate for constructing 
a stable MaVaN scenario. 
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I. INTRODUCTION 



Substantial observational evidence arising from Type-la supernova data [l|-l3|, big-bang 
nucleosynthesis constraints j4|, cosmic microwave background radiation (CMBR) power 
spectrum ^7|, large scale structure js], and determinations of the matter density jo], Q 
suggest a cosmological model where the energy density (p) of the universe is dominated by 
an overall smoothly distributed component with negative pressure (p), the so-called dark 
energy, which leads to an accelerated expansion. 

At the current stage of our knowledge, speculations about the nature of the dark energy 
and the negative pressure component responsible for the accelerated expansion of the uni- 
verse have given origin to a fruitful discussion in the literature llMlSl]. The most obvious 
explanation for dark energy is the cosmological constant, a constant vacuum energy density 
which has as equation of state, p = —p, at all times. However, since the cosmological con- 
stant has a magnitude completely different from that predicted by theoretical considerations, 
physicists have been compelled to consider other explanations for the dark sector 16l-l22|. 
A plausible alternative for obtaining a negative pressure equation of state, better motivated 



by the high energy physics, considers a dynamica 
field rolling down in a fairly flat potential 



23 



energy density governed by a light scalar 



25| . These models assume that the vacuum 



energy can vary [26]. Guided by theoretical as well as phenomenological arguments, several 



possibilities have been proposed, such as /c-essence 
and also several types of modifications of gravity 



31 



33|. 



phantom energy models 29l. |30|. 



In this context, one of the most challenging issues concerns models of mass varying 
neutrinos (MaVaN's) 34-3 

coupled to the dark energy light scalar field component. The 
simplest realization of the MaVaN mechanism consists in writing down an effective potential 
which, in addition to a scalar field dependent term, contains a term related to the neutrino 
energy density. 

On the other hand, dependin g on an eventual compatibility with well-known cosmological 



bounds on neutrino masses 



38 



411 ] ■ for some choices of neutrino-scalar field couplings and 



scalar field potential, which is intermediated by the so-called stationary condition 36|], the 
combined fluid (dark energy plus neutrinos) is subject to instabi 



become non-relativistic (NR). Due to the instability problem 37|, |42|, models so far proposed 



ities once the neutrinos 



usually face flne-tuning problems, which impose severe constraints on the choice of the 
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scalar field potential and their kinetic energy 43|]. Indeed, when a scalar field potential is 
tuned in agreement with the abovementioned stationary condition, an effective potential, 
which includes a cosmological term dynamically dependent on the scalar field and a neutrino 
contribution, is supposed to have a minimum with a steep second derivative for a finite scalar 
field vacuum expectation value 37|]. If this condition is not satisfied, an attractive force 
between neutrinos intermediated by the scalar field can lead to the formation of neutrino 

n 

nuggets [44]. An immediate consequence is that the coupled fiuid would behave as cold dark 
matter and not suitable dark energy candidate. 

Notice that dark matter has not been considered in the formulation of the MaVaN mech- 
anisms. However, the idea of unifying dark energy and dark matter naturally offers this 
possibility. The generalized Chaplygin gas (GCG) [iJQ is particularly relevant in this 
respect as it is shown to be consistent with the observational constraints from CMB 45| . 



supernova 
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48| . gravitational lensing surveys j49|, and gamma ray bursts [50]. Moreover, 



it has been shown that the GCG model can be accommodated within the standard structure 
formation mechanism 
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In this paper we consider the possibility of neutrino masses arising from an interac- 
tion with a real scalar field which describes the dynamics of the generalized Chaplygin gas 



171 . |47| ]. As is well known, the Higgs sector |5l| and the neutrino sector are possibly the 



only ones where one can couple a new standard model (SM) singlet without upsetting the 
known phenomenology. As it shall be seen, one can treat the mass varying neutrino term 
in the equation of energy conservation as a perturbative contribution to the evolution of 
the previously unperturbed adiabatic solution, as for obvious phenomenological reasons, we 
expect a small contribution from neutrinos to the energy dynamics of the universe. In par- 
ticular, we suggest that the dynamics of the scalar field is modified by a linear perturbation 
over the field itself due to the neutrino mass coupling, which is actually turned on when the 
neutrinos become NR. 

The basic steps for introducing this perturbative approach as well as the procedure for its 
implementation is presented in section II. The applicability of this procedure is, on general 
terms, compared with the applicability of the stationary condition in the extended Fardon- 
Nelson-Weiner (FNW) framework |36|. In section III, we review the main features of the 
GCG, its underlying scalar field (0), and its dependence on the universe scale factor, a. 
In section IV, we discuss the neutrino mass generation mechanism in the context of the 
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GCG model. For this purpose we consider left-handed non-sterile neutrino masses with an 
analytical dependence on 1/0 and 1/0^ as suggested by the simplest version of the see-saw 
mechanism. In section V, we set the adequate phenomenological constraints for the neutrino 

masses and for the GCG parameters, consistent with the most recent observational bounds. 
Finally, in section VI, we draw our conclusions. 



II. A PERTURBATIVE APPROACH FOR COSMOLOGICAL NEUTRINOS 

Gosmological neutrinos have not yet been observed, so hints about their nature and 
behaviour are quite welcome. In the usual MaVaN framework, neutrinos are coupled to a 
light scalar field which is identified with the dark sector. Presumably, the neutrino mass rrii, 
has its origin on the vacuum expectation value (VEV) of the scalar field and, naturally, its 
behaviour is governed by the dependence of the scalar field on the scale factor. Since the 
neutrino statistical distribution corresponds to a Fermi-Dirac statistics without a chemical 
potential f{q), where q = Ti,^ being the neutrino background temperature at present, 
the neutrino energy density and pressure can be expressed in the following way 

PMc^) = 4^ Cdq^' U + ^^^y^ (1) 



where we have assumed a flat FRW cosmology and introduced the sub-index for denoting 
present-day values, with ao = 1. For simplicity, wc have considered a single non- vanishing 
neutrino mass, although the generalization to more than one massive neutrino is straight- 
forward. 

Simple mathematical manipulation allows one to easily demonstrate that 

m„{(j)) ^"^"''f^^ ^ {pu{a, (/)) - Sp^{a, (/>)), (2) 

where, for the neutrino NR regime, 

— ~ n^{a) (xa . 3) 

From the dependence of on a, one can obtain the energy-momentum conservation for the 
neutrino fluid, 

p^{a, cj,) + 3H{p^{a, <j)) + p^{a, </>)) = , (4) 



where H = a/a is the expansion rate of the universe and the overdot denotes differentiation 
with respect to time (' = d/dt). 

It is important to emphasize that the couphng between cosmological neutrinos and the 

■ i)q. (|21) is restricted to times when neutrinos are NR, i. e. 



scalar field as described by 
(Ty(a) = T^o/a » m^{(j){a))) , the decoupled fluids should evolve adiabatically since the 



nj,(a) oc a 



36 



371 . |43| . In opposition, as long as neutrinos are relativistic 



strength of the coupling is suppressed by the relativistic increase of pressure (p^ ~ Sp^). In 
this case, one would have 

Pu,,p + SH{p^^^ + p^^^) = 0, (5) 
where for the scalar field background fluid. 



(6) 



Treating the system of NR neutrinos and the scalar field as a unified fluid (UF) adia- 
batically expanding with energy density pi,p = + Pcf, and pressure PuF = Pu + P<f) lead 
to 

puF + 3H{puF + Puf) = P0 + 3H{p^ + p^) = -4>—^t;z^, (7) 



d(f) drriy 

where the last step is derived from the substitution of Eq. (jl]) into the first equation in (JTj). 

The existence of a scalar field dark energy sector on its own constitutes a problem in what 
concerns order of magnitude equality with the energy densities of the other components of 
the universe. The theoretical assumptions proposed in Ref. Jsoj] and subsequently developed 



by other authors 



22 



37 



42 



43| suggest a stationary condition which allows circumventing 



the coincidence problem for cosmological neutrinos in a way that the dark energy is always 
diluted at the same rate as the neutrino fluid, that is, 



dV{(t>) duly dpy 
dd) d(h duly 



(8) 



This introduces a constraint over the neutrino mass since it promotes it into a dynamical 



quantity, as indicated in Eq. (JTj). The main feature of this scenario [36| is that it is equivalent 
as adopting cosmological constant like dark sector, with an energy density that varies as a 
function of the neutrino mass. As already mentioned, the effectiveness of this coupling is 
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restricted to values of the scale factor larger than a^vB, where af^u parameterizes the transition 
between the relativistic and NR regimes. 

The assumption of a universe with the dark energy sector governed by the equation of 
state p0 = — p0 implies, through Eq. (|H]), that p\ = V, and allows to recover the stationary 



condition, exactly as obtained in Ref. [36|. It is clear that the relevance and the con- 
siderations about the constraints on the equation of state are actually dominated by the 
competition among scalar field potentials and its adequacy for neutrino mass generation. 
In fact, once one assumes that = — p^, the neutrino mass evolution and the form of the 
potential become automatically entangled by the stationary condition. Thus, it should be 
reahzed that the stationary constraint Eq. ([S]) is quite dependent on the potential of the 
scalar field. 

Moreover, the introduction of a kinetic energy component modifies the equation of state 
and implies that the stationary condition is not satisfied when Pd) + = 4>^ is non- vanishing. 



This difficulty has already been pointed out in Ref. [43||, where a solution which severally 
constrains the choice of the scalar field potential is proposed. It is then shown that the FNW 
scenario is consistent only for a vanishing kinetic energy contribution, for a dark energy fiuid 
behaving like a "running cosmological constant". 

It follows from this discussion that, at least from a theoretical point of view, an alterna- 
tive approach to treat deviations from the FNW proposal is needed to overcome the quite 
restrictive condition ptf, + P4, = 0, which, in turn, implies in a vanishing kinetic energy con- 
tribution. In spite of being, at present, negligible in comparison to the potential terms, this 
contribution can dominate and affect the evolution of the universe at earlier times. Scalar 
field dark sector candidates with a well defined equation of state, such as for instance the 
GCG, are natural alternatives to circumvent this problem. Candidate models can be easily 
implemented by means of an assumption that the neutrino coupling to the underlying scalar 
field of the dark energy sector is perturbative. In this instance, we start instead of Eq. (I?!), 
from the unperturbed equation (JSj), and establish the conditions for treating the neutrino 
coupling in a perturbative way. To exemplify this point, let us consider the unperturbed 
equation of motion for the scalar field, 

and assume that the effect of the coupling of the neutrino fluid to the scalar fleld fluid is 



6 



quantified by a linear perturbation e0 (|e| << 1) such that 



(j) ^ (p {1 + e)(j). (10) 

It then follows the novel equation for the energy conservation 

y I 3gy I ^^^^^ - (11) 
dtp dip dniy 

The explicit dependence of on is easy to quantify. Indeed, after a simple manipulation 



one finds 



dV{ip) 
dip 



_^dV{^) 



dW{cl>) 



dV{^) 



(12) 



The substitution of the Eqs. (ITOl) and (fT2l) into Eq. (fTTj) and use of Eq. (|9]), lead to 



(13) 



where the perturbative character of the neutrino mass term is assumed when we set the last 
approximation in the above equation. Finally, we can obtain for the value of the coefficient 
of the perturbation 

(14) 



d(j> (9m 1/ 



d ( 1 



which for consistency is required to satisfy the condition |e| << 1. Of course, this procedure 
can be carried out to higher order as to determine the additional perturbative modifications 

to PuF- 

It is important to notice that in the conventional quintessence models the scalar field 
is, at the present epoch, slowly rolling dawn its potential and therefore its effective mass, 
[dW / d(l)'^y^'^ is smaller than the Hubble expansion rate, H. In opposition, the usual SC 
framework for MaVaN's establishes that the dynamical scalar field sits at the instantaneous 
minimum of its potential, and the cosmic expansion only modulates this minimum through 
changes in the local neutrino density [4^. It hence allows for {d^V / dcjy^Y^^ » H, which 
means that the coherence length of the dynamical scalar field is much smaller than the 



present Hubble length, and thus, unlike in quintessence models, the perturbations on sub- 
Hubble scales are adiabatic. Consequently, the speed of sound, for these perturbations is 
simply given by cl = dp^p/dp^p. The very argument of a scalar potential being modulated 
by a local neutrino energy density applies to our approach. This is reinforced by the idea 
that the Higgs fsi] and the neutrino sectors are the only ones where one can couple a new 
standard model (SM) singlet without upsetting the known phenomenology. 

Otherwise, our setup allows one to consider a large class of scalar field potentials and 
equations of state for the dark sector, for which the stationary condition is incompatible 
with generic neutrino mass generation models. Furthermore, given its perturbative nature, 
our procedure is valid for any equation of state satisfying the relevant phenomenological re- 
quirements. In the following analysis one can verify the stability condition by observing that 
the square of the speed of sound of the coupled fluid is also dominated by the unperturbed 
equation of state. 

In fact, this feature and the ensued problem could have already been observed in the 
context of the FNW scenario since the stationary condition implies that = — 1 and the 
burden of recovering the stability (c^^^^^^ > 0) is transferred to the neutrino contribution. 

It is easy to see that our result reduces to the FNW scenario when = = or when 
one assumes that = —p^ = V. Indeed, rewriting the equation for the conservation of 
energy, Eq. ([7]), by simply redefining p^p as 



the usual definition 
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PUF = -0' + Kff, (15) 

33 for an effective potential Kfi in terms of ^ = 



is now valid for any value of (p. 

In the next section we shall implement this perturbative approach and verify the applica- 
bility criteria for the coupling between variable mass neutrinos and the scalar field associated 
with the GCG background fluid. 



III. THE GENERALIZED CHAPLYGIN GAS - AN STATE EQUATION FOR THE 
DARK SECTOR 



The GCG model is characterized by an exotic equation of state [19|, |45| given by. 



p = -A,po ( — 
P 



(16) 



which can be obtained from a generahzed Born-Infeld action 19|. In any case, irrespective of 
its origin, several studies yield convincing evidence that the GCG scenario is a phenomeno- 
logically consistent approach to explain the accelerated expansion of the universe. The 
constants and a are positive and < a < 1. Of course, a = corresponds to the ACDM 
model and we are assumin g th at the GCG model has an underlying scalar field, actually 



real 



0,0 



or complex 



18, loi. The case a 



1 corresponds to t 



le equation of state of the 



Chaplygin gas scenario 17] and is already ruled out by data j45|. Notice that for As = 0, 
GCG behaves always as matter whereas for A^ = 1, it behaves always as a cosmological 
constant. Hence to use it as a unified candidate for dark matter and dark energy one has to 
exclude these two possibilities so that Ag must lie in the range < A^ < 1. 

Inserting the above equation of state into the unperturbed energy conservation Eq. (E]), 



one obtains through a straightforward integration 19 1 

:i-Ar 



P<f> = Po 



A, + 



3(i+a) 



i/(i+a) 



and 



P<p = -A,po 



A, + 



fl- A 



Following Ref. |47[, one can obtain through Eq. (j6]) the field time-dependence. 



Po(l - A, 



a 



3(a + l) 



As + 



A„ 



3(i+q) 



a 



and assuming a flat evolving universe described by the Friedmann equation 
H in units of and p<^ in units of pcrit = SH^JSnG), one obtains 



0(a) 

where one assumes that 



2/3 



l_A,(l-a^/3) + v/r^ 



00 = 0(ao = 1) 



1 



In 



1- 



i + yr^j 



(17) 



(18) 



(19) 
(with 



(20) 



(21) 



and P = 3(a + l)/2. 

One then readily find the scalar field potential in terms of the field, 

V{cj>) = ^Aj^po {[cosh (/30)]^ + [cosh (;g0)]"^| . (22) 

The analytical dependence of the energy density, the potential energy and the scalar field 
in terms of the scale factor is illustrated in the Fig. [1] for some particular choices of a and 



for = 0.7. The latter value arises from the matching of the model with supernova, CMB 



data, and cosmic topology 



45 



47 



52|. Notice that in what concerns a, the observational 
constraints are as follows: WMAPl is compatible with a < 0.6 45|; WMAP3 admits values 



in the range a < 0.2 53|]; and structure formation implies that a < 0.2 

One of the most striking features of the GCG fluid is that its energy density interpolates 
between a dust dominated phase, pch oc a~^, in the past, and a de-Sitter phase, pch = —Pch, at 
late times. This property makes the GCG model an interesting candidate for the unification 
of dark matter and dark energy. Indeed, it can be shown that the GCG model admits 
inhomogeneities and that, in particular, in the context of the Zeldovich approximation, 
these evolve in a qualitatively similar way as in the ACDM model 19|. Furthermore, this 
evolution is controlled by the model parameters, a and A^. 

In order to understand the possible range of values for a, one has to consider the prop- 
agation of sound through the GCG fluid, = dp^/dp^. A detailed quantitative analysis of 
the stability conditions for the GCG in terms of the squared speed of sound is discussed in 
Ref. [47|. Positive cl implies that < a < 1. At a later stage, we shall come back to this 
issue when considering the coupling with neutrinos. 



IV. NEUTRINO MASS MODELS 



Despite its impressive phenomenological success, it is widely believed that the SM of par- 
ticle physics is actually only a low-energy approximation of an underlying more fundamental 
theory. In this respect, the interplay with the cosmology can be an important guideline to 
obtain insights on the nature of the more fundamental theory. In the SM, the most natural 
way to explain the smallness of the neutrino masses is through the seesaw mechanism, ac- 
cording to which, the tiny masses, m,^, of the usual left-handed neutrinos are obtained via 
a very massive, M, sterile right-handed neutrino. The Lagrangian density that describes 
the simplest version of the seesaw mechanism through the Yukawa coupling between a light 
scalar field and a single neutrino flavour is given by 

C = mLRULUji + M{(t))uRUR + h.c, (23) 
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where it is shown that at scales well below the right-handed neutrino mass, one has the 



effective Lagrangian density 



55| 



Phenomenological consistency with the SM implies that logarithm corrections to the above 
terms are small, while it is well-know from the results of solar, atmospheric, reactor and 
accelerator neutrino oscillation experiments that neutrino masses given by 171,^ ~ m\^^/ M{4>) 



lie in the sub-eV^ range 56| • It is also clear that promoting the scalar field into a dynamical 



quantity, (j) — )■ ip{a) leads to a mechanism in the context of which neutrino masses are time- 
dependent. Associating the scalar field to the dark energy field allows linking NR neutrino 



energy densities to late cosmological times [22|, I36|,l37^. This scenario can be implemented 



through the perturbative approach via Eq. ([7]). It is evident that this approach is fairly 
general, as well as independent of the choice of the equation of state and on the dependence 
of the neutrino mass on the scalar field. For sure, the form of M{4>) and of the equation 
of state can lead to quite different scenarios. In studying the neutrino coupling with the 
GCG fluid, we consider two cases, namely, rriy = mo{(j)a/(t>) and rrii, = mo{(f)a/(t>y. These 
two scenarios are by no means the only possibilities. In particular, we choose them as they 
correspond to the simplest feasible possibility for sterile neutrino effective mass generation, 



55|. 



given that the scalar field has mass dimension one 

In particular, the essential information content of the mass dependence on the scale factor 
is obvious from the fact that neutrino masses increases with decreasing 0, which in the GCG 
model is a decreasing function of a. In the phenomenological analysis performed in the next 
section we shall elucidate this point by noticing that for larger values of the scale factor 
(decreasing red-shift) the neutrino mass increases. 

In addition to the mass dependence, it is necessary to determine for which values of the 
scale factor the neutrino-scalar field coupling becomes important. For convenience we set 
the value of a = apfn for which p,y^MR = Pu,ur (which sometimes is expressed in terms of 
rriu ~ T^), to parameterize the transition between the NR and the ultra-relativistic (UR) 
regime. In fact, this occurs when 

m,ia) = mo(0o/0(a))" = X— (25) 

where n = 1, 2 and x is a numerical factor estimated to be about x ~ 94 supposing that 
Pu/Pciit = rrio [ey]/(94 /i^ [eV]), where h is the value of the Hubble constant in terms of 

11 



100 km s~^Mpc~^ . This condition allows one to establish the correspondence between the 
values of a^vij for which the (UR to NR) transition takes place and the neutrino masses 
assume the present-day values. For our model such a transition is shown in Fig. |2] for the 
studied set of parameters of the GCG model. 

The consistency of our perturbative approach depends on the maximal value assumed by 
the e parameter at present (cf. Eq. (fl^ ). In order to estimate e, we use Fig. [31 From Eq. ( 120|) 
we observe that the dominant analytical dependencies of (p and V{(j)) on the scale factor can, 
in good approximation, be expressed as 0(a) ~ K^\n{a) and V{4){a)) ^ + K-^a'^^"^^"^ in 
the interval < a < 1, where arbitrary constants. These lead to a 

dependence of e on the scale factor given approximately by 

°^ ln"+^(a) (ln(a) + K,)' ^^^^ 
This is an increasing function of the scale factor dominated by the mass dependent term 
contribution, m,^(0(a)) oc 0(a)~", n = 1, 2. By noticing that eMax = e(a = 1), we show in Fig. H] 
the dependence of eMax for which GCG parameters and ttIq values are consistent with our 
perturbative approach. 

In the next section we analyze the phenomenologically interesting set of parameters for 
our GCG-neutrino coupled model. As mentioned, the choice of = 0.7 will be considered 
as it is consistent with all known data. 



V. THE STABILITY CONDITION AS A WINDOW FOR THE NEUTRINO PHE- 
NOMENOLOGY 



Current cosmological data constrain the number of active neutrino flavours as well as 

3 

the sum of their masses to m^j < 0.75 eV at 95% c.l. [38i]. This constraint does not 

i=i 

agree with the Heidelberg-Moscow bounds arising from of neutrinoless double beta decay 



which sets the limit ^^m^j > 1.2 eV at 95% c.l. 57|, which, however, are in agreement 

i=i 

with the CMB data analysis of the WMAP results that sets 0.7 eV at 95% c. 1. for each 
neutrino species (and 2.0 eV in total) 58j. Actually, improvement on experimental data are 



expected to be sensitive to the effects of a finite sum of neutrino masses as small as 0.06 eV 



38 



59| . the lower limit arising by neutrino oscillation experiments that set Am^ ~ 7.5 



8.7 X 10 eV^ (2a) for solar neutrinos, and Am^ ~ 1.7 - 2.9 x 10 ^ eV^ {2a) for atmospheric 
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neutrinos. In addition, current neutrino direct mass measurements, for instance, through 
the tritium beta decay Mainz experiment, does set an upper hmit on the effective electron 
neutrino mass of < 2.3 eV at 95% c.L [60]. Thus, in order to test our scenario, we shall 
consider present neutrino masses varying from 0.05 to 5eV. 

the GCG parameters using supernova and CMB data as 



On the other hand, studies o 



well as cosmic topology 



45 



47 



471 . |52| allow choosing a typical value = 0.7, for about 



70% of dark energy and 25% of dark matter in a universe filled with 95% by the GCG fluid. 

With these values we can illustrate the behaviour of the neutrino mass in terms of the 
GCG parameters as a function the scale factor. This is depicted in Fig. [5l One sees that 
for ~ 1/0 neutrinos become NR earlier than the rrij, ~ 1/0^ case. Actually, the stronger 
is the inverse power dependence of the neutrino mass on the scalar field, that is m,^ ~ (f)'", 
n > 1, the later neutrinos become NR. If on the other hand, the mass generation model were 
exponential type, mj, ~ exp— [0/0o], then the neutrinos would become NR much closer to 
present. 

In addition, from the constraint set by Eq. ( I25l) . which establishes the neutrino NR 
regime, we show in Fig. El for the same set of parameters, present-day values of the neutrino 
masses and the corresponding values of ajvjj for which the transition between the NR and UR 
regimes takes place. Since there are strong phenomenological constraints on the choice of 
TTLo, it is important to pay attention to the present-day value of neutrino mass interval, from 
0.05 to 5eV, where a clear dependence on the model for mass generation is observed. 
The clnr values found in Fig. |6]are taken into account in Fig. [TJ where it is examined the 
validity of the perturbative approach for different GCG scenarios. 

By observing the model dependent conditions discussed above, namely the value of a = 
ajvij, and the present-day neutrino mass corresponding to the maximal value of the linear 
perturbation coefficient eMax, we can, for instance, set the phenomenologically acceptable 
values of mo = 1 eV and mo = 0.1 eV, in order to perturbatively quantify the modifications 
to the energy density components of the coupled fluid. Such perturbative corrections are 
shown in Fig. [HI 

Interestingly, for mo = 1 eV, a fairly typical value, we can see in Fig. [9] that stable 
MaVaN perturbations correspond, for the GCG case, to a well defined effective squared 
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speed of sound speed 

^ > 0. (27) 

The larger are the mo values, the larger are the corrections to the squared speed of sound, 
up to the limit where the perturbative approach cannot be applied. Therefore, as far as 
the perturbative approach is valid, our model does not run into stability problems in the 
NR neutrino regime. In opposition, in the usual treatment where neutrinos are just coupled 
to dark energy, cosmic expansion together with the gravitational drag due to cold dark 
matter have a major impact on the stability of MaVaN models. Usually, for a general fluid 
for which we know the equation of state, the dominant effect on c\ arises from the dark 
sector component and not by the neutrino component. For the models where the stationary 
condition (cf. Eq. ([H])) imphes a cosmological constant type equation of state, = —pip, 
one obtains = — 1 from the very start of the analysis. For sure the situation cannot be 
fixed by the perturbative contribution of neutrinos. Our GCG-neutrino model is free from 
this inconsistency. 

VI. CONCLUSIONS 

In this work we have analyzed the coupling of neutrino masses to a light scalar field 
associated with the GCG unification model of dark energy and dark matter. We have 
presented the criteria of applicability of a perturbative approach for in the study of MaVaN 
models and determined the coefficient of a linear perturbation which is given in terms of 
the ratio between the variation of the neutrino energy and scalar field potential terms (cf. 
Eq. f|T4l) ). Seesaw proposals for the neutrino mass were then considered and the effects of the 
neutrino mass coupling with the underlying scalar field of the GCG model were quantified. 
As discussed, our approach yields a positive squared speed of sound and is consistent with 
the current neutrino mass experimental limits. 

For sure, adopting a perturbative approach is equivalent to the assertion that the coupling 
between neutrinos and dark sector is fairly weak. Besides, our proposal turns out to be an 



interesting alternative to the usual stationary condition constraint proposed in the Ref. [36 1 
for the equation of state = —p^. For several configurations, the latter scenario leads to 
catastrophic instabilities associated with an imaginary speed of sound in the neutrino NR 
regime. Actually, in previous work it has been pointed out that MaVaN models generically 
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face stability problems for some choices of neutrino-scalar field couplings and scalar field 
potentials once neutrinos become NR. Effectively, the scalar field mediates an attractive 
force between neutrinos which can lead to the formation of neutrino nuggets. This would 
turn the combined fluid to behave like cold dark matter and thus render it non-viable as a 
candidate for dark energy. 

In opposition, our analysis shows that the coupling of neutrinos with the scalar field of 
the GCG model is consistent in what concerns the positiveness of the squared speed of sound 
(cf. Fig. [9]). The knowledge of the background fluid equation of state and our perturbative 
approach allows one to overcome the problematic negative squared speed of sound. Even 
though our analysis could be extended to other equations of state, the GCG model seems 
to be a quite suitable candidate for constructing stable MaVaN scenarios. 

To conclude, we stress that allowing for dynamical behaviour to a scalar field associated 
to dark energy in connection with the SM neutrinos and the electroweak interactions may 
bring important insights on the physics beyond the SM. Neutrino cosmology, in particular, 
is a fascinating example where salient questions concerning SM particle phenomenology can 
be addressed and hopefully better understood. We believe that our proposal is a further 
step in this respect. 
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FIG. 1: The energy density p<^(a), the potential energy V{(j)) and the scalar field (/>(a) as functions 
of the scale factor a for several choices of a and for A, = 0.7. 
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FIG. 2: Present values of the neutrino mass mo and the corresponding values of a^R for which 
the transition between the NR and UR regimes takes place in the GCG phenomenological sce- 
nario with variable As and a = 1, 1/2. The neutrino mass varies with 1/0 and l/cp'^ and the 
increasing gray level corresponds to increasing values of mo, for which the boundary values are 
mo = 0.05 eF, 0.1 eV, 0.5 eV, leV, 5eV. 
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FIG. 3: Determination of the maximal value of the |e| parameter as a function of the scale factor 
and of the GCG parameters and a. For illustration have set mo = 1 eV . The increasing 
gray level corresponds to growing values of |e|, for which we have marked the boundary values 
|e| = 1, 0.1, 0.01. Independently of the neutrino mass dependence on (j) (m or m ~ it 

is clear from the graphs that eMax = e(a = 1) for all that cases we have considered. 
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FIG. 4: Maximal value of the linear perturbation coefficient e as a function of the present-day 
neutrino mass and of the energy density coefficient for the GCG. Once again we have set the 
neutrino mass varying as l/(^ and 1/0^ in a GCG scenario for a = 1, 1/2. The increasing gray 
level corresponds to increasing values of mo, for which we have marked the boundary values for 
mo = 0.05 eF, 0.1 eV, 0.5 eV, leV, 5eV. 
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FIG. 5: The neutrino mass {mjj{a)jm^ dependence on the scale factor a. We have set model 
dependent neutrino masses varying as 1/^ and Xj^"^ in a GCG scenario with = 0.7 and a = 
1, 1/2, 1/3, 10-^ 
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FIG. 6: Present-day values of the neutrino mass m-o and the corresponding values of Ojvii for which 
the transition between the NR and UR regimes takes place in a GCG phenomenological scenario 
with = 0.7 and a = 1, 1/2, 1/3, 10"". The choice of the model for mass generation plays a 
relevant role in determining the starting point ajva of the coupling effectiveness. This is a section 
of graphs of Fig. [2] for = 0.7. 
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FIG. 7: The present-day neutrino mass rrto with respect to the maximal value of the linear per- 
turbation coefficient e for = 0.7 in a GCG scenario with a = 1, 1/2, 1/3, 10"*. As before, the 
choice of the model for mass generation plays a relevant role in determining the range of validity 
of the perturbative approach. Notice that the maximal value for e corresponds to a = 1. 
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FIG. 8: Components of the energy density (in units of pcmt) obtained from our perturbative 
approach for = 0.7 in a GCG scenario with a = 1, 1/2, 1/3, 10"". The choice of the model for 
mass generation does not play a relevant role on behaviour of the modified energy components. 
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FIG. 9: Perturbative modification on the square of the speed of sound c^(a) as a function of the 
scale factor for the neutrino-GCG coupled fluid in comparison with the adiabatic GCG fluid for 
As = 0.7. The couphng with the GCG turn the MaVaN's naturally stable. 
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